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SUMMARY 


The response to atmospheric gusts of a lightly damped airplane 
model flying at a Mach number of 0.7 at sea level has been studied analyt- 
ically with elevator fixed and with varying amounts of viscous restraint 
of a mass -overbalanced elevator. Although in a sharp-edge gust the vis- 
cously restrained mass -overbalanced elevator has negligible effect on the 
motion of the airplane center of gravity until the first peak in the 
normal-acceleration response of the model has been reached, the effective 
damping ratio of the subsequent motion can be more than trebled by a suit- 
able choice of viscous restraint with some reduction in frequency of 
oscillation. In continuous gusts, calculations indicate that for the 
model considered there results a 20-percent reduction in root-mean-square 
normal acceleration at the center of gravity over a broad range of eleva- 
tor viscous restraint. 


INTRODUCTION 


Aircraft operating in gusty air experience loads which can be divided 
into two types, according to their source: (l) loads due to the direct 

effect of the gust and (2) loads due to motions induced by the direct 
gust loads. Loads of the first type increase with airspeed, and air den- 
sity (see ref. 1), and those of the second type are significantly depend- 
ent on the damping ratio of the short -period oscillation of the airplane 
(see ref. 2) and the airplane relative -density factor (see ref. 1). 

With lightly loaded well-damped airplanes, the motions induced by 
the gusts serve to reduce considerably the overall loads In rough air; 
however, with heavily loaded lightly damped airplanes, the reduction due 
to the motions is much less. Thus, means of gust-loads reduction become 
increasingly more desirable as the airplane relative-density factor 
Increases and the short-period damping decreases. This reduction can be 
achieved (l) by countering the direct effect of the gust (see, for exam- 
ple, refs. 1 and 3 to 10) and (2) for airplanes with low damping ratios. 
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by altering the motions induced by the direct gust loads by increasing 
the short-period damping ratio. To date there has been relatively little 
work in this second direction, perhaps because of the limited alleviation 
that might be ejqpected from such a system (see ref. 11). 

The purpose of the present paper is to investigate analytically the 
alleviation that can be obtained by a relatively simple system which 
alters the longitudinal motions induced by the direct effect of the gust 
loads. This simple system consists of a mass-overbalanced viscously 
restrained elevator and is analogous to the simplified yaw damper reported 
in reference 12. The calculations presented in this paper are limited to 
an airplane model flying at one altitude and at one speed. This airplane 
model has a relatively large pitching radius of gyration and static margin. 
Calculations were made for a wide range of viscous restraint, and the 
resulting reductions in normal acceleration at the center of gravity were 
determined. 


ANALYSIS 


Description of System and Outline of Analysis 

The system studied in this paper is an airplane with a mass- 
overbalanced viscously restrained elevator and is similar to the simpli- 
fied yaw-damper system reported in reference 12, except for the use of 
mass overbalance in lieu of adjusted hinge -moment parameters. It was 
decided to use elevator mass overbalance in this investigation because 
of its inherent simplicity and the difficulty in predicting the effect 
of elevator modifications on hinge moment. The system has seven degrees 
of freedom for the rigid airplane, but when the motions are restricted to 
the longitudinal modes with essentially constant forward speed, only 
three degrees of freedom remain: vertical translation of the center of 

gravity of the airplane, rotation of the airplane about the lateral axis 
through the center of gravity, and rotation of the elevator about its 
hinge axis. In the following analysis, the differential equations 
governing the motions of this system in its three degrees of freedom are 
set up, and the dynamic stability of the system is studied from the roots 
of the characteristic equation. Then, the response of the system to a 
unit step gust input (indicial response function) is obtained and used to 
obtain the response to continuous sinusoidal turbulence (the frequency- 
response function) . Finally, the root-mean-square normal acceleration 
in continuous random turbulence is obtained by utilizing power- spectral- 
density techniques. 
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Equations of Motion 

In the analysis the usual assumptions pertinent to the longitudinal 
dynamic stability of a rigid aircraft with free elevators were made. The 
resulting differential equations of motion, based on the body-axis system 
shown in figure 1, are 

Fg = mV(<x - 9) = qS (pz a a + (l) 


M y = IyG = qSc (c^a, + f + | C^.fl + C^S + | C^.b) (2) 


Mft = I e (9 + 6) + l e m e x e 0 - nigXgVCe - a) 

= <lS e 5 e (c^a. + ^ C^a + ^ C^.8 + 0^8 + C^.s) (3) 

(The symbols are defined in appendix A.) The force Fg is related 
to through the equation 

•»--$ w 

The characteristic equation obtained from these differential equa- 
tions is given in appendix B and has the form 

?\(aA^ + BA 5 + CX 2 + DA + e) =0 (5) 

The characteristic equation (5) is of fifth degree with one zero 
root, and there are three distinct types of motions, depending on the 
value of c hg,* If the value of is less than a certain critical 

value, there result two sets of complex conjugate roots corresponding 
to the oscillatory airframe and elevator modes and the one zero root 
which represents the indifference of the airplane to pitch attitude. 

If c hg, exceeds the critical value, two unequal negative real roots 

result, representing the exponential decay of the elevator mode; the 
airframe mode is still represented by the complex conjugate pair and 
the zero root. If equals the critical value, the two negative 
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roots are equal. This condition corresponds to critical damping of the 
elevator mode. This special case is not dealt with in this paper. 

i/ 


Response to Unit Step Gust Input 

An elaborate formulation of the problem of the response of the system 
to a unit step gust input, such as that given by Mazelsky and Diederich 
(ref. 13) and Brenner and Isakson (ref. 14) is not given herein. Instead, 
the normal-acceleration indicial-response model proposed by Preps and 
Mazelsky (ref. 2) was used in a modified form. In this normal-acceleration 
response model it is assumed that after the airplane encounters a step gust 
directed normal to the flight path “the shape of the response curve up to 
the peak acceleration depends primarily on the gust penetration function 
(Kttssner function) . After the peak value . . . the character of the response 
appears to be primarily a function of airplane stability and to approximate 
the short-period oscillation of the airplane. 11 (See ref. 2.) Or, in equa- 
tion form. 


a n (s) = a njmax sin (0 < s < s p ) (6) 


Ctq_ ( s “Sp) 

®nC 8 ) = a n,max e cos K s ' s p) (s P < s > 


This response model has been used successfully by Vitale, Press, and 
Shufflebarger (ref. 15) in predicting root-mean-square normal acceleration 
on a transonic tailless rocket -propelled model. In the present study a 
more nearly exact representation of the short-period oscillation was 
used and the penetration effect on the tail was accounted for in an 
approximate manner that agrees at least qualitatively with the tail effects 
shown by Mazelsky and Diederich (ref. 13) and Brenner and Isakson (ref. 14) . 
In reference 14, it is seen that the tail begins to develop an appreciable 
fraction of its quasi-steady lift only after it penetrates several tail 
chords into the gust and that, for high mass ratios, practically no transla- 
tion or pitching occurs until this time. Accordingly, it was assumed that 

(1) No vertical translation or pitching of the airplane occurred 
until after the tail had penetrated 2 tail chords into the gust (total 
penetration of 4 wing chords). 

(2) The wing lift equaled its quasi-steady value after it had pene- 
trated 4 wing chords into the gust and was a sine function of the pene- 
tration distance before this time. 
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(3) The tail lift equaled its quasi-steady value after it had pene- 
trated 2 tail chords into the gust and was zero before this time. 

The assumptions are not valid for all airplanes, but are thought to 
apply sufficiently well for this configuration in a heavily loaded condi- 
tion (high mass ratio) . The discontinuous development of the tail lift 
causes discontinuities in normal acceleration and pitching moment which 
approximate the more exact variations indicated by reference 14 for an 
airplane with a high mass ratio. 

The elevator effects were treated by assuming that 

(1) The transient aerodynamic hinge moments occurring up to the first 
peak were insignificant in comparison with the artificial viscous hinge 
moments and the inertia moments of the mass -overbalanced elevator. 

(2) Subsequent to the first peak, the aerodynamic hinge moments 
could be represented adequately by quasi-steady values. 

The normal-acceleration and pitching responses up to 4 wing chords 
penetration can be written down immediately as a consequence of the 

assumptions. Thus, since Og = 


F Z _ 

TT - 


/-C 


Za,WB 


■i) sin 


m n 


( 8 ) 


0=0 


(9) 


The elevator motion obtained by solving equations (l) and ( 3 ) 
simultaneously subject to the conditions of equations (8) and (9) and 
of zero aerodynamic hinge moment but with artificial viscous restraint 
is given by 


( 10 ) 


As was pointed out previously, it is assumed that, subsequent to the 
first acceleration peak, the airplane is free of unsteady- lift effects 
due to penetration into the gust; that is, its lift can be given ade- 
quately by quasi-steady values. This assumption is equivalent to the 
assumption that the model responds essentially in its short-period longi- 
tudinal motion subject to the initial conditions 


6 = Kq 


-1 


coJC, 
L 0 1 



UJ. &t , , c-1 

—2. e x - sin cd t + — = cos 


ov 
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a ( t p) " 7 

8(t p ) = 0 

S(tp) = 6p , 

S(t p ) = 0 

8(t p ) - 5p 


( 11 ) 


where Bp and Bp are obtained from equation (10) at t = tp. The solu- 
tions of equations (l), (2), and (3) subject to the initial conditions 
given by equations (ll) are the solutions required. These solutions, 
representing an for the two types of motions, are given as functions of 
the nondimens ional penetration distance s = 

c 


( 12 ) 


For c h 6 » < ( C h & .) cr -’ 

r 

a n( s ) = ij - jA x cos p(s - s p ) + B ± sin p(s - s p )J + 

^ ■> 
e a 2( p) ^ cqs A ^ s _ ^ + sin a(s - B p ^l 

J J 

For Ch 5' > ( Ch 6 ') cr’ 

a n( s ) = fj |e ai ^ B " S p) Ja x cos p (s - s p ) + sin p (s - s p )J 

C ie a 3 ( s -s p ) + Die o4(s-s p )J 


( 13 ) 


For , = -oo (elevator fixed), 


a n( s ) = ‘ e ai ( S_Bp )[A 1 cos p (s - s p ) + B 1 sin p (s - s p ) 


( 14 ) 


V 




u 


to 
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Hie expressions for 0 and 6 are similar and are given in appendix C 
together with the procedure for the determination of the constants A-^, 

B;j_, C]_, and Dj_. 

Note that, for the elevator-fixed case, the representation of the 
short-period mode (eq. (l4) ) contains a sine term which does not appear 
in the simplified representation assumed by Press and Mazelsky (eq.. (7))- 
Although the resulting representation of the short-period mode (eq. (l4)), 
is more nearly exact, it is only an approximate representation of the 
oscillatory response of the airplane subsequent to the first peak, since 
the penetration effects on the aerodynamic derivatives probably persist 
for more than the 4 chords assumed. 


Response to Continuous Gusts 

The gust frequency-response function or the response of the airplane 
to an infinitely long succession of gusts of the form 


a g = i sin (B g t 


(15) 


can be obtained from the normal-acceleration indicial-response functions 
(eqs. (12) and ( 13 )) by the relation 


r 

F n (iv) = iv / a n (s)e“ 1 

Jo 


-rvs 


ds 


(16) 


(Xw-C 

where v = The integrations for the two cases have been performed 

and are given in appendix D. For subsequent calculations with the atmos- 

owc 

pheric power-spectral-density function, the frequency variable v = 
is transformed to 0 = ^£. Thus n is independent of the airplane speed 
and size and is based solely on atmospheric characteristics. 

The power spectrum of the normal acceleration was ccaiputed from the 
equation 


4 n< n > - 4 at m os< a >K< ln > 


(17) 


The power-spectral-density function of the atmosphere was assumed to be 
given by 
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p 

(intensity constant) i /,p\ 

4 atmo S ( Q ) 52 - 5? (I8) 

In the present study, variations in atmospheric gust intensity were not 
considered and the unit intensity constant was used for all calculations. 

The root -mean- square normal acceleration was then found from the 
relation 


ff n = 



ll/2 


$ n (n)dn 


(19) 


These results (eqs. (17) to (19)) from the generalized harmonic analysis 
of gust loads are taken from reference 2. 


CALCULATIONS AND RESULTS 


An outline drawing of the airplane configuration chosen for applica- 
tion of the foregoing theory is given in figure 2. The straight horizontal 
tail was chosen to give a straight hinge line. Some of the pertinent char- 
acteristics of this airplane are given in table I. The size was chosen to 
be comparable to the rocket -powered models tested by Vitale, Press, and 
Shufflebarger (ref. 13 )• The aerodynamic characteristics of the wing-body 
combination were based on tests by Spreemann and Alford (ref. l6)j the tail 
and elevator data were taken from results of a test by Turner (ref. 17)* 

When this investigation was first undertaken, the roots of the charac- 
teristic equation were obtained for this airplane with a mass-balanced 
elevator; the results indicated relatively little effect of the elevator 
mode on the airframe mode for any value of C^g , • A value of mass over- 
balance (xg negative) sufficient to halve the E term in the characteristic 
equation was chosen for further calculations . For lg flight under the 
assumed flight conditions, the trim elevator angle corresponding to the 
mass overbalance was -0.47°. 

The period and damping of both the airframe and the elevator short- 
period modes are given in figure 3 as functions of C^g, . (Note that 

C^g, = -oo corresponds to elevator fixed. ) The gust indicial responses 

in a n and 6 are shown in figure 4 for Chgt = 0* -62.5 > -175 > and -°°j 

the 8, § , and S responses are also given for Qhg, = -175 and -oo. The 

remaining calculations were done for Cj^, = 0, -62.5, -175 > -250 > and -oo. 
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The amplitude ratios of the normal-acceleration gust-frequency-response 
functions and the corresponding power-spectral-density functions are shown 
in figures 5 snd 6, respectively. The root -mean- square normal-acceleration 
response (fig. 7) was then determined by graphically integrating $ n (fl). 


DISCUSSION 


Stick-Free Stability 

The effect of viscous elevator restraint on the classical stick-free 
stability of the airplane is shown in figure 3- It is seen that the addi- 
tion of viscous damping can more than treble the effective damping of the 
airplane short- period mode (as compared with stick-fixed), and is extremely 
effective in improving the damping of the elevator mode as could be expected. 
(Note that Ch§i due to aerodynamic effects alone is only on the order of 

-1 to -10.) In the vicinity of Qhg, = -80, the aforementioned critical 
value of Cjjg, is reached. At this point the osci llat ory elevator mode 

disappears and is replaced by two exponential convergences; for further 
increases in Chg, the motion can be practically represented by a single 

exponential decay. It is interesting to observe that the peak in the 
airframe -period curve occurs in the vicinity of ( C hgi) cr - 


Gust Response 

Discrete gusts .- Comparison of the time histories of the gust indicial 
responses for Ch 5 , = 0 , - 62 . 5 , -175 , and -«> (fig. 4) shows that viscously 

restraining the elevator markedly influences the motion subsequent to the 
first peak. For Chg, = -175 > £ a = O. 565 , and the subsequent oscillations 

practically disappear in one cycle; whereas for Chg t = -«, £ a = O. 175 , and 

about one-third of the original amplitude remains after 1 cycle. This 
favorable effect of the elevator motion is accomplished with only small 
elevator deflections (6 = 0 . 05 °) and velocities (6 = 2 . 2 ° per second) for 
a gust velocity of 1 foot per second. Hence the drag penalty due to the 
action of the elevators would be insignificant. 

Continuous gusts .- The normal-acceleration gust-frequency-response 
function for Qhg t =0 (fig- 5 (a)) exhibits two peaks; the low-frequency 

peak corresponds to the oscillatory airframe mode and the higher frequency 
peak to the oscillatory elevator mode. In considering the magnitude of 
the elevator-mode peak, it must be kept ■ in mind that there is no damping 
at all in the elevator circuit, and the addition. of even a slight amount 
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of viscous restraint corresponding to the aerodynamic damping of the 
elevator (ph§ t = ’ bo " 10 ) would - substantially lower this peak (see 

fig. 5(b)). The small amount of damping which exists in the elevator v 

mode is provided by the dynamic coupling with the airframe mode. For 
Cftg , = -62.5, although the elevator motion is still oscillatory, the peak 

has disappeared. For values of v > 0.2, the frequency-response function 
is practically the same for all values of Cgg t because of the predomi- 
nating contribution of the initial quarter sine wave in the normal- 
acceleration indicial response at these frequencies. 

The power spectrum of the normal-acceleration response (fig. 6) indi- 
cates that practically all the power in the normal-acceleration response 
occurs near the natural frequency of the airframe. Note that, for 

= -175, the normal acceleration at the airframe natural frequency 

has been reduced 48 percent below that experienced with Cgg, = - 00 but 
that at lower frequencies the acceleration is somewhat greater. 

The principal results of this investigation are presented in fig- 
ure 7. Therein, the root-mean-square normal acceleration of the airplane 
in continuous rough air, normalized with respect to that obtained with 
fixed elevator, is given for four values of vIbcous restraint. The * 

ratio a n ^b n 5 fixed is seen to be about 1.07 for no viscous restraint 

and decreases rapidly with the initial increase in Cg & , • A flat minimum 

of about 20-percent alleviation is reached for Cg^, = -125 to -250; 

although computations were not made for higher values of Cg^, , it seems 

reasonable to expect that the curve would follow a gradual upward trend 
with increasing Cg^, and approach 1.0 asymptotically. 

The dashed line in figure 7 labeled "Simple indicial- response model" 
is based on the simplified response model used by Press and Mazelsky 
(eqs. (6) and (7)) and permits the determination of cr n if the effective * 

values of airplane short-period stability parameters, that is, damping 
ratio an< * f re Q. uenc y b a , be known. The omission of the airframe sine 

term from the indicial response model (see eqs. (12) to (l4)) evidently 
causes the simplified method to overestimate the alleviation near optim um 
Cggi and to underestimate it for low values. 

Other considerations .- In all the foregoing analysis and application, 
the moments acting on the elevator have been assumed to be significantly 
larger than moments due to static friction. Also, no provision was made 
for elevator control, since the elevator used for the analysis occupied 
the entire exposed trailing edge of the tail. In a practical application 
to an actual airplane, it might be feasible to employ a split elevator, 


v 
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one-half being used for control and one-half being used for auxiliary- 
damping. Such an arrangement would also insure that the elevator motions 
required for damping would not be interfered with by the action of the 
pilot in holding the control stick. The system could be optimumized at 
cruising speeds; if the action of the system proved objectionable at the 
lower speeds experienced during climb and descent, it could be locked 
out, since the highest speeds likely to be attained are the critical 
speeds . 

The effects of changes in speed, Mach n umb er, and altitude on the 
optimum values of Phgt have not been investigated and the optimum values 

of I e , Xg, Cj^, and for this and other values of static margin 

have not been determined. Also, since mass overbalance is used, there 
would be a small hinge moment due to longitudinal accelerations and a 
reduction in stability in steady pull-ups, with possibly sane effect on 
the phugoid mode. Xn spite of these possible effects, the system does 
appear to offer a usable margin of gust -load alleviation for lightly 
damped aircraft, and it is especially attractive in view of its inherent 
simplicity and reliability. 


CONCLUDING REMARKS 


The gust response of a lightly damped airplane model flying at a 
Mach number of 0.7 at sea level has been studied analytically with eleva- 
tor fixed and with varying amounts of viscous restraint of the mass- 
overbalanced elevator. Although the viscously restrained mass-overbalanced 
elevator probably has negligible effect on the airplane motion until the 
first peak in the normal-acceleration response to a sharp-edge gust has 
been reached, the effective damping ratio of the subsequent motion can be 
more than trebled by a suitable choice of viscous restraint with some 
reduction in frequency of oscillation. In continuous gusts, calculations 
indicate that for the model considered there results a 20-percent reduc- 
tion in root-mean-square normal acceleration at the center of gravity over 
a broad range of elevator viscous restraint. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., September 10, 1957* 
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APPENDIX A 
SYMBOLS 


Positive directions of displacements, velocities, and accelerations 
are shown in figure 1. 


a real part of characteristic -equation root pertaining to 

airframe mode 

normal -acceleration increment from trim due to gust of 1 ft/sec, 
g units 

b imaginary part of characteristic -equation root pertaining to 

airframe mode 

c real part of characteristic -equation root pertaining to eleva- 

tor root 

£ mean aerodynamic chord of total wing 

c e mean aerodynamic chord of exposed elevator 


, - ^ 
h qS e c e 



F 7 

C 7 = -g 
L qS 

d imaginary part of characteristic-equation root pertaining to 

elevator mode 

e base of natural logarithms, 2. 7182 . . . 

F n (iv) ,F n (ift) gust frequency-response functions 
Fg force parallel to Z-axis 

g acceleration due to gravity, 32.17 ft/sec^ 
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i = \pi 





L e,2 


k. , - SsSbl 

* J-e “i 


mass moment of inertia of elevator . about hinge axis 


c , -Jk. 

qSc, 


e,2 ~e,l 5, 


I- , i 

c e 


mass moment of inertia of airplane about lateral axis through 
center of gravity 




= Jx 

qSc 


r t V 

E Y,2 * £ 


Z e distance from center of gravity of model to center of gravity 

of elevator 


m 


mass of airplane 


m, =S§ 

1 qS 


% 

My 


mass of elevator 

aerodynamic hinge moment on elevator taken about hinge axis 
pitching moment about center of gravity 
free-stream Mach number 


P 


period of oscillatory motion, sec 
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<1 

dynamic pressure, lb/sq, ft 

i* 

6 

distance penetrated into gust, taken from leading edge of 
exposed wing root chord, 

P 

s 

wing area 


t 

time 


\/2 

time for motion to damp to one -half amplitude 


u 

component of V along X-axis 


V 

free-stream airspeed 


w 

weight 


w 

component of V along Z-axis 


x e 

distance from elevator hinge axis to center of gravity of 
elevator, positive when center of gravity is rearward of 
hinge axis 

V 

X 

longitudinal body axis, positive forward 


Y 

lateral axis, positive toward right wing tip 


Z 

normal body axis, positive downward 


a 

angle of attack, radians 


„ _ ac 

°i "T 



H 

^loi 



A*c 
“5" V 



«, -M 

4 V 



a g = ? 


H 



V 
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o _ be 

P “ *7“ 


7 

5 


a dc 

A = T 


C 
e 

V = 


03gC 

V 


atmos 




(n) 


CD, 


8 


flight -path angle , radians 

elevator displacement, radians (unless otherwise noted) 


damping ratio 
pitch attitude, deg 

root -mean-square normal acceleration 

power- spectral- density function of the atmosphere 

power-spectral-density function of the normal- acceleration 
response in continuous rough air 

circular frequency of sinusoidal gusts 


CD a = 


2s t 


CD = 

° 2t T 


a -3a 

v 

Subscripts : 


a 

cr 

e 

WB 

P 


airframe 
critical 
elevator 
wing -body 

pertaining to the first peak in the normal-acceleration response 


^fixed 


elevator fixed 
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Derivatives are denoted by subscript notation, thusly: 


dC 

C<7 = — — and so forth 

^ da 


Differentiation vith respect to time is 


denoted by a dot, thusly: 


a 


da 

dt 


whereas differentiation with respect to distances (in wing chord lengths) 
is denoted by a prime, thusly j 


a * = £££ = a _ da c 
ds v/S dt‘ V 


however , 


9 


, _ d& c 

dt T 


K 


b 
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APPENDIX B 


THE CHARACTERISTIC EQUATION 

The characteristic equation obtained from the equations of motion 


is 


a(aA^ + BA 5 + CA 2 + DA + e) = 0 


(Bl) 


where 


A = 1 


(B2) 


B 


„ ~ Ch 5' _ f Cm Q' + + £Za\ + L + z e m e x e\ Ss 


L e,2 


V-Y,2 I Y J 2 m l 


I e ] X Y,2 


(B3) 


= + + . %_ + £ + S) . % + 

I ej2\lY,2 %,2 m l j X Y,1 \ X e / X Y,1 X e,l 


% 

CniQ t 

L . z e®eXe^ 

| Cln B' 

+ 

/l + l e m e x e\ 

m l 

X Y,2 

V Ie ) 

1 J Y,2 

+ m l 

LV *•) 


L Y,2 


" G hs' _ Vm^Xg\ 

X e y 


°niB * rka 1 . * 

T 


^,2^,2 Z e,2> 


(B4) 
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D = 


C Za Sa\ 

4- 

K 

“l "%,!/ 

T _* 

“l 

X e,l 


' Ch 5 A f° h 8» ^ VmeXe^^i 


< J e t 2 I e J J X,2 


f 1 + 1 e m e x e\ Cm 5 




r \k* _ Vffl e*e\ S) » 
^e,2 I e 7^,2 


m i 


( l e mgxA 0 ^ c h a /C h VfogX^St' 

V 1 "T77 + — J + 

+ + SxA _ ^°6 ( V + \ _ Ch a 

I e,l\ I Y,2 ^2/ I Y,1\ I «,2 ^,2/ X Y,2 I e ,l 


(B5) 


E = 


/Gina, Chg 

Cms Cha \ 

♦ 5i[ 

/Ch 0 , ( Vm e x e > 

C ®6 

i 

(O 

£ 

1 

1 

VY,1 X e,l 

X *A I e,l/ 

“x 

ye, 2 I e j 


X Y,2 x e,lJ 


CZ; 


'5 


mi 


C hg C me' _ / C h0 '_ + Vm e x e\ St 


^,1 *1,2 l x e,2 ' J e J 1 Y,1 


(B6) 


The roots of the characteristic equation are as follows: 
For Chj., < (Chg,)^, 

*1,2 = a ± ib 

which corresponds to the short-period mode of the airframe j 


Aj ^ = c ± id 

which corresponds to the short-period mode of the elevator^ arid 


Aip = 0 

which corresponds to indifference of the airframe to pitch attitude. 
For C hs , > (Ch 5 ,) cr , 
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A-j ^ “ a i ib 

which corresponds to the short -period mode of the airframe j Aj and 

are negative real roots which correspond to the exponential decay of the 
elevator mode; and 

A^ = 0 

which corresponds to pitch indifference. 
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APPENDIX C 


RESPONSE TO UNIT STEP GUST INPUT 


The responses to a unit step gust input are given in terms of the 
nondimens ional penetration distance s = 


For 0 < s < s_, 


-c z , 


ajs) = - ■ 1 sin £ 

ct 2 a~ 


g 


(ci) 


0 ( s ) = 0 


S(s) = 


-1 


03 o K- 


o*l 


1 




+ K^AKl 


Kqc 
o e — 


0> 0 c Ko 
sin ■ ■ s + —i 1 
V B * <% 


0> o C 

-ri COS -X- 


(C2) 


(C3) 


For s p < s < co and C hg . f < 


r 

( “Cv \ a i( s “ s p) r 1 

“5^ « e [ A l cos ^ ( s - Bp) + B 1 sln p( s - Sp)! 

<x 2( s ~ s p)r r \-i 

e Cp cos A^s - s p ^ + D-l sin A(s - s p J 


(C4) 


I ccq_ f s — s p ^ p 

0(s) = \e A2 cos p(s - Sp) + B£ sin p (s - s p ) + 


a 2 ^s ” s p ) p \ t “j 

i C 2 cos A^s - s p ) + D 2 sin A^s - s p ) 




+ e 2 (C5) 
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5 (s) = e a l( s-s p) [a 3 cos p(s - s p J + sin p(s - Spj| 

a 2( s “ s p)r i \ t \\ 

e C5 cos A(s - SpJ + sin Afs - SpJ 

For s p < s < co and C h5I > (Ch 5 .) cr * 


(06) 


M 8 ) = 


r -C, 


Z ^Ne 


Oi(s-Sp) 


\ m i e 

0l e"5( a - s p) + Dl e^( S - S P)| 


[aj cos p fs - SpJ + sin p(s - Sp)l 


( 07 ) 


a l( s-s p) r 

d(s) = e A 2 cos P ^s - Spj + B 2 sin P (s - SpH + 

a, (s-s ) cc 4 (s-6p) 

C 2 e 5 5 + D 2 e ' + E 2 


(08) 


Op (s-Sp) |- 

8 (s) = e |a 5 cos p(s - Sp^ + B^ sin p (s - s p ) + 


“3 ( S_S P) ^ t-> 
5 T 


C^e ■*" 


( 09 ) 


It should be pointed out that the arbitrary constants Eire A-j_, Bi, C]_, 
I> 2 _, and E 2 . The remaining constants, A 2 , A^, B 2 , B^, Cg, C^, D 2 , 

and are the modal constants and are linear functions of the arbitrary 

constants and the stability derivatives, mass parameters, and flight con- 
ditions. These constants were determined by standard mathematical 
techniques . 
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APPENDIX D 


RESPONSE TO CONTINUOUS GUSTS 


The normal-acceleration response of the airplane to an infinitely 
long succession of gusts of the form 


1 1 /vV\ 

cc,g = ^ sin “g'k - y sinl-E-Jt 

is given by the following equation for C^g, < (Ch§r) cr : 


(Dl) 


B'n(iv) CZq, WB l 
m lS 


“s' 


2 -gjv cos vs p + l(u E - v Bin vSp)]| 


\v^ + 2(0^2 - p 2 )y2 + (aj 2 + p 2 ) 


jvv - p 2 + v 2 j cos vsp + (-<* 1 ) + P 2 + v 2 ) sin vsj 

|" v ^ x l 2 “ P 2 + v 2 jsin vs p + (-Oj^c^ 2 + P 2 + v 2 ) cos VBpj^ + 

V* * a(^ - f' “* * K 2 * p2 - ->] * 

jj-v^-2ax)sin vs p + ^aj 2 + p 2 - V 2 ) cob vSpj^ + 

C-ij it ■■ vfao 2 - A 2 + v 2 )cos vs_ + { -a^A (a., 2 + A 2 + v 2 )sin vs_l 

Y + 2 (^ 2 - + 1? - *) 1 L ^ p 1 pJ 

i^ 2 - A 2 + v 2 ) sin vsp + ^-ag^ctg 2 + A 2 + v 2 jcos vsjj 

d i(t — n — ifr~ 7 2 — iT5([ v (- 2a 2) c ° 8 vs p + .(“ 2 2 + a2 - v2 ) sin "»] + 

\v 4 + 2 fag 2 - A 2 J v 2 + (ag + A 2 y 
i |jv(-2a 2 ^sin vSp + (a 2 2 + A 2 - v 2 )cos V3pJ^j 


(D2) 
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For c hgi > ( c hgi) cr > the expression for F n (iv) is modified only in the 
Cp and terms. These terms now become 

g • [ v cos vs p + (-o^sin vspj + i jj-v sin vs p + (-a^j cos V8 p]^ 


2 ][ V C ° S VS P + (-^ii) 13111 vs p] + i Q-v sin vs p + (-cc^)cos vSp" ► 

\ ^ ^ Jj 


In either case, the initial term 


2 2 

V - v 


Jv cos vs p + i^cn s - v sin vSp)J 


= JL + i 
4 2 


when cd — v. 

o 
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TABLE I.- PERTINENT AIRPLANE CHARACTERISTICS 


Geometric characteristics: 

S, sq ft 

c, ft 

Aspect ratio 

Sweepback, deg 

Taper ratio, deg 

Distance from center of 
gravity to 0.25 mean 
aerodynamic chord, ft . . . 


Wing 

Stabilizer 

Elevator® 

3.77 

O.565 

0.10 

1.068 

0.584 

0.0895 

4.0 

4.0 


b 45.0 " 

b 7.13 

c 0 

0.500 

0.600 

O.665 

0 

2.51 

d 2.50 


Mass characteristics: 


m, slugs 5.00 

nig, slugs 0.0264 

Iy, Slug-ft^ 6.60 

Iy e , slug-ft 2 0.000295 

*e/ c e °*5366 


Flight conditions: 

Altitude, ft 0 

Static pressure, lb/sq ft 2,ll6 

Velocity of sound, ft /sec 1,116 

0-7 

q^, lb/sq ft 725 

V, ft /sec 780 


Aerodynamic characteristics: 


c z a 

. -4.25 



. . . -4.15 

% 

. 0 

Omg 1 ........ 

. . . 0 

Cz a,WB 

• 

1 

VJJ 

. 

Co 

-3 

Ch a 

. . . -0.100 

Sp, wing chords 

• 4 

ChB 

. . . -0.484 

Sx 

-1.162 

Ch a' 

. . . -1.95 



-0.556 

C 4e* 

. . . -5.28 

%x* 

-0.87 

C 4b< • 

. . 0 to -a> 

a 

All calcinations based 

on exposed elevator area. 



13 At quarter -chord line. ■ 

c 

At hinge line. 

distance from center of gravity to elevator hinge axis. 




Figure 1. - Body-axis system used in analysis. Axes remain fixed relative to airplane. 
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Figure 4.- Continued. 
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(c) C^, = -175. 


Figure 5.- Continued. 
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Figure 6.- Continued 
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Figure 7-~ Normalized root-mean-square value of normal-acceleration response as affected by- 

variations in elevator viscous restraint. 


--5 



